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We consider the charged-current quasielastic scattering of muon neutrinos on an Oxygen 16 target, described within a 
relativistic shell model and, for comparison, the relativistic Fermi gas. Final state interactions are described in the distorted 
wave impulse approximation, using both a relativistic mean field potential and a relativistic optical potential, with and without 
imaginary part. We present results for inclusive cross sections at fixed neutrino energies in the range E v = 200 MeV - 1 GeV, 
showing that final state interaction effects can remain sizable even at large energies. 



1. INTRODUCTION AND FORMALISM 

A precise and realistic description of neutrino- 
nucleus interaction in the intermediate energy region 
is crucial for the interpretation of experimental results 
used to determine neutrino properties. At present 
most of the Monte Carlo codes 1 1 1 which have been 
developed to simulate the response of the detectors in 
these experiments are based on the Fermi gas model. 
Therefore they take into account the Fermi motion of 
the nucleons inside the nucleus and Pauli blocking ef- 
fects, but they neglect several other effects, which are 
known to be important from electron scattering exper- 
iments [ 2 1 . 

In order to study the relevance of some of these 
effects, in this contribution we compare relativistic 
Fermi gas (RFG) calculations of neutrino-nucleus 
quasielastic scattering cross sections J3J, with results 
obtained within a relativistic shell model (RSM), in- 
cluding, in particular, final state interactions (FSI) [4|. 
More details about our calculations can be found in 
Ref. |5|. FSI effects in v-nucleus scattering have also 
been studied in Refs. 1617 1 . 

We consider the charged-current (CC) quasielastic 
scattering of muon neutrinos on 16 O at fixed neutrino 
energies of 200, 500 and 1000 MeV, describing this 
process within the impulse approximation, schemati- 
cally represented in fig.[D 

We thus assume that the incident neutrino interacts 



with only one nucleon, which is then emitted, while 
the remaining (A-l) nucleons in the target are spec- 
tators, that the nuclear current is the sum of single 
nucleon currents and that the states of the target and 
residual nuclei are described by independent particle 
model wave functions. We describe the ground state 
of 16 as a closed shell configuration, the occupied 
shells being s^, P3/2 an d Pi/2- F° r the removal of 
a nucleon from a closed shell of angular momentum 
j, the cross section corresponding to the diagram in 
fig Eh as the following general form: 

x %^El"( fe ')Y a (l+Y 5 )«W-/a(q)| 2 ^A-i, 

where Gf is the Fermi constant, £ indicates the aver- 
age/sum over the initial/final spins, 75 = — /^y'^y 3 
and Dirac spinors are normalized according to 
u(k)'u(k) = 2Icq. We then sum over the occupied 
shells and integrate over the emitted nucleon and over 
the direction of the outgoing muon in order to get the 
inclusive cross sections do/dT^, and the integrated in- 
clusive cross section o = J (do / dT M ) dT^, T p being 

the outgoing lepton kinetic energy. 

The main ingredient in eq. Q is the single nucleon 
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Figure 1. Born approximation diagram for CC V- 
nucleus quasielastic scattering. The impulse approxi- 
mation is assumed in the hadronic vertex and the pos- 
sibility of final state interactions between the outgoing 
nucleon and the residual nuclear system is explicitly 
shown. 



current matrix element, 

-/a(q) = VvJ rfV' qr y^,r)f„<'(r) , (2) 

where \\tg\r) and V|/. SA ,(pw, r) are the wave functions 
for the initial (bound) nucleon and for the emitted nu- 
cleon, respectively, and F a is the single nucleon weak 
charged-current operator. For the latter we assume the 
free, on mass shell, nucleon expression: 
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where Fyja are the CC single nucleon Pauli and 
Dirac form factors, Fa and Fp are the axial and the 
induced pseudoscalar form factor, respectively, and 
\V M i\ is the ud Cabibbo-Kobayashi-Maskawa matrix 
element. For Fa we assume a dipole parameterization 
with cutoff mass M A = 1 .026 GeV. 

In eq. (0 the bound nucleon wave functions (r) 
are the self-consistent (Hartree) solutions of a Dirac 
equation, derived, within a relativistic mean field ap- 
proach, from a Lagrangian containing o, co and p 
mesons, which has been already successfully used in 
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Figure 2. Differential cross section (da/dT M ) versus 
the outgoing muon kinetic energy, for the quasielastic 
scattering of muon neutrinos on 16 O and incident neu- 
trino energy: E v = 200 MeV (upper panel), 500 MeV 
(middle) and 1 GeV (lower panel). 



the study of (e,e'p) processes |4|. For the single- 
particle binding energies of the different shells we use 
the corresponding experimental values, which deter- 
mine the threshold of the cross section for every shell. 

Concerning the emitted nucleon, as a starting point 
we describe it as a plane wave, thus neglecting its in- 
teraction with the residual nucleus (plane wave im- 
pulse approximation - PWIA). Then, to make our de- 
scription more realistic, we include the effects of FSI 
by using distorted waves (distorted wave impulse ap- 
proximation - DWIA). For the latter, we make the fol- 
lowing different choices. 

Complex ROP: following previous studies of ex- 
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Figure 3. Same as fig.|2j but including FSI effects. All 
curves are calculated using the RSM model, in PWIA 
(solid) and describing FSI within the RMF (dashed), 
real ROP (dotted) and complex ROP (dot-dashed) ap- 
proaches. 



elusive electron scattering processes |4], we employ 
distorted waves which are obtained as solutions of a 
Dirac equation containing a phenomenological rela- 
tivistic optical potential (ROP). The ROP has a real 
part, which describes the rescattering of the ejected 
nucleon and an imaginary part, that accounts for the 
absorption of it into unobserved channels. 

Real ROP: since, contrary to the (e,e'p) case, here 
we are considering inclusive processes, where all fi- 
nal channels contribute, the presence of the imaginary 
term in the optical potential leads to an overestimation 
of FSI effects. For this reason we also consider the 
potential obtained by setting the imaginary part of the 



ROP to zero (a discussion on the use of real optical 
potentials has been presented in Ref. |7 1). 

RMF: finally, we employ wave functions which are 
obtained as the solutions in the continuum of the same 
Dirac equation which is used to derive the bound nu- 
cleon wave functions. We refer to this approach as 
relativistic mean field (RMF) and consider it appro- 
priate at low energy transfer. 

2. RESULTS AND CONCLUSIONS 

Let us now present the results we obtain for 
neutrino-nucleus quasielastic cross sections. As a first 
step we neglect FSI effects and just compare RSM- 
PWIA results with the corresponding curves obtained 
within the relativistic Fermi gas, for which we use a 
Fermi momentum pp = 225 MeV and binding energy 
eg = or 20 MeV. This comparison is presented in 
fig- El where the differential cross sections da/dT^ is 
plotted as a function of the outgoing muon kinetic en- 
ergy. We observe that the differences between the two 
models are quite large at low neutrino energy, but they 
practically disappear at E v — 1 GeV. 

The situation is different if we "turn on" FSI ef- 
fects, as illustrated in fig. where the RSM-PWIA 
results of fig. [2] are compared with the DWIA ap- 
proaches previously outlined. 

We see that FSI effects produce a reduction of the 
cross section, with respect to PWIA. The RMF and 
real ROP curves are quite similar, showing a reduction 
of about 30 4- 40% for E v = 200 MeV and 20% for the 
other energy values. In the case of the complex ROP 
model, instead, the reduction of the cross section is 
rather large, between 60% (E v = 200 MeV) and 50% 
(at higher energy), due to the absorption introduced 
by the imaginary term. 

Let us note that the resonant structure appearing at 
high Tfj is due to the use of real potentials (RMF and 
real ROP) for describing the final nucleon state. Here 
we only include single-particle excitations within a 
mean field picture, while including residual interac- 
tions would make the width and number of resonances 
to be considerably larger. 

Finally, the results of figures |2] and [3] are summa- 
rized in fig. |4] where the cross section o, integrated 
over the muon energy, is plotted as a function of 
the incident neutrino energy E v . Again we see that 
within the PWIA the discrepancy between different 
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Figure 4. Integrated cross section a for the quasi- 
elastic scattering of muon neutrinos on 16 as a func- 
tion of the incident neutrino energy. The curves are 
calculated within the RFG model, while the points 
correspond to RSM calculations without FSI (stars) 
and with FSI effects taken into account within the 
RMF (empty squares), real ROP (full squares) and 
complex ROP (circles) approaches. 



nuclear models is relatively small and decreases with 
increasing neutrino energy, while FSI effects are still 
present even at large E v . The (more reliable) results 
for the RMF and real ROP approaches show a not too 
large, but still appreciable reduction (~ 15% at E v = 1 
GeV), while, again, the imaginary term in the com- 
plex ROP leads to a too large reduction (~ 50%) of 
the cross section. 

In conclusion, within the Impulse Approximation, 
we observe that when FSI effects are neglected, 
the nuclear model dependence of the quasielastic v- 
nucleus cross section is large at low neutrino energy 
but becomes negligible at E v = 1 GeV. FSI effects are 
also rather large at E v ~ 200 MeV and decrease with 
increasing neutrino energy, but, in the more realistic 
real ROP approach, can still be as large as 15 % at E v 
= 1 GeV and thus must be carefully considered. 

As a final remark, we notice that if we consider 
Neutral Current quasielastic processes, where the out- 
going nucleon must be detected, then using the com- 
plex ROP description of FSI would be more appropri- 
ate. Therefore in this case rather large FSI effects are 
expected to be present even at relatively high neutrino 



energy I3I8I . 
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